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In quantum chemistry, the price paid by all known efficient model chemistries is either the truncation of the Hilbert space or 
uncontrolled approximations. Theoretical computer science suggests that these restrictions are not mere shortcomings of the 
algorithm designers and programmers but could stem from the inherent difficulty of simulating quantum systems. Extensions of 
y—{ computer science and information processing exploiting quantum mechanics has led to new ways of understanding the ultimate 
<0> limitations of computational power. Interestingly, this perspective helps us understand widely used model chemistries in a new 
light. In this article, the fundamentals of computational complexity will be reviewed and motivated from the vantage point of 
(3JT) chemistry. Then recent results from the computational complexity literature regarding common model chemistries including 
^ Hartree-Fock and density functional theory are discussed. 
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\£> 1 Introduction 

Quantum chemistry is often concerned with solving the 
' ' Schrodinger equation for chemically-relevant systems such as 
atoms, molecules, or nanoparticles. By solving a differential 
■ and/or eigenvalue equation, the properties of the system and 
Q the dynamics of the state are obtained. Examples of proper- 
ties include: equilibrium geometries, the dissociation energy 
of molecules, and the vibrational frequencies. 

The difficulty stems from the accuracy required and the ap- 
q parent exponential growth of the computational cost with both 
" ^ the number of electrons and the quality of the description of 
^ the system. For practical applications, the accuracy required is 
typically orders of magnitude smaller than the total energy of 
the system. As a concrete example, the Carbon atom has total 
electronic energy of about 37.8 Hartrees while the energy of 
t— I a Carbon-Hydrogen bond is only 0.16 Hartrees. Solving the 
full eigenvalue equation takes on the order of n 3 operations 
for an n x n matrix. However, when describing interacting 
many-electron systems, the dimension of the matrix increases 
exponentially with the number of electrons. 

Consequently, the computational methods of electronic 
structure in chemistry are aimed at circumventing exact di- 
^vq agonalization, in the context of electron structure, called the 
t-H full configuration method. Avoiding exact diagonalization has 
^ led to a wide range of computational methods^B] f or com- 
• <"H puting properties of chemical interest. These methods have 
recently^] begun to include quantum simulation following 
^ Feynman's suggestion^ to use quantum computers as sim- 
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ulators. Subsequent development of these ideas in quantum 
chemistry has led to new proposed methods utilizing quantum 
computational technique s^ 11 " 14 ! and proof -of -principle exper- 
iments 15 16 . However, questions about when and w here one 
would expect a quantum computer to be useful^E! have not 
been fully answered. Ideally, computational complexity can 
provide some answers about when, where, and why quantum 
computers would be useful for chemistry. At the same time, 
it could also help formalize intuitive understanding of when 
we can expect reliable results from traditional computational 
methods. 

Many results in computational complexity have inter- 
changed classical and quantum computers fluidly leading to 
new results on the complexity of computing properties of 
quantum systems 19 . We review some recent results appearing 
in the computational complexity literature that touch on why 
electronic structure calculations are difficult. Our hope is to 
encourage future investigations into quantitative understand- 
ing of difficult instances of electronic structure calculations. 

A similar, but much shorter, discussion of computa- 
tional complexity in quantum chemistry by Rassolov and 
Garashchuk 20 appeared in 2008 and provided many conjec- 
tures that have since been proven or extended. 

This perspective assumes exposure to second quantization 
and mixed states in quantum mechanics. Standard bra-ket no- 
tation is used when referring to quantum states. All the nec- 
essary concepts from computational complexity and quantum 
computation are briefly introduced and motivated to make the 
article as self-contained as possible. A key omission from the 
pres ent re view is a number of computational complexity re- 
sults EE3 concerning quantum-information based wave func- 
tion ansatzes such matrix product states, density matrix renor- 
malization group and their generalizations 23 . These methods 
are becoming accepted into the mainstream of computational 
chemistry 24 " 26 ! b u t do not yet have the widespread availability 
of the selected methods included in the present work. 
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2 Worst case computational complexity for 
chemists 

The purpose of this section is to provide a quick but precise 
introduction to computational complexity conceptst 2 ^ 2 ^. The 
aim is to set the stage for the results subsequently reviewed. 

Computational complexity is the study of how resources re- 
quired to solve a problem change with its size. For instance, 
space complexity is the scaling of memory requirements with 
the problem size, but, in this article, we focus on time com- 
plexity which investigates how the running time of the compu- 
tation changes as the problem size increases. 

Computational chemists often informally discuss common- 
place concepts in computer science such as the com- 
plexity classes of polynomial-time problems (P) and non- 
deterministic polynomial-time problems (NP). For instance, 
it is sometimes stated that Hartree-Fock has a runtime which 
scales as the third power of the number of basis functions. 
This is scaling disregards difficult instances of the calculation 
where Hartree-Fock does not converge. Such instances require 
manual intervention to tweak the algorithm used or adjust the 
convergence thresholds in a case-by-case fashion. 

In computer science, often the focus is on worst-case com- 
plexity where the most difficult instances of a problem are used 
to classify the problem's complexity. This is has been one 
of the major areas of theoretical computer s cienc e, although 
work on average case complexity does exist 29 30 . In this ar- 
ticle, the complexity of the problems discussed are character- 
ized by the worst-case scaling. 

2.1 Time complexity in equivalent computer models 

A proper measure of the time complexity of an algorithm is 
how many basic operations (or how much time) it takes to 
solve problems of increasing size. Conventionally, a compu- 
tational problem is described as easy or tractable if there exists 
an efficient algorithm for solving it, i.e. one that scales poly- 
nomially, 0(n k ) with fixed k and input size n|^| Otherwise, the 
problem is considered intractable. This is an asymptotic def- 
inition that may not capture the full utility of the algorithm. 
For example, an asymptotically efficient algorithm may run 
slower than an asymptotically inefficient for small or fixed 
size problem instances. Nevertheless, this asymptotic clas- 
sification of algorithms has proved useful. From a theoretical 
computer science perspective, the division allows for consid- 
erable progress to be made without considering the minutiae 
of the specific system, implementation, or domain of applica- 
tion. 



* The notation, f(x) = 0(g(x)), implies that function f[x) is dominated by g(x) 
for asymptotically large values of x. To indicated that fix) is asymptotically 
larger than g(x), we write f(x) = Q.(g(x)). If f(x) is dominated by and dom- 
inates function g(x), then we write f[x) = &(g(x)). 



From a practical standpoint, Moore's^ law states the den- 
sity of transistors in classical computers doubles every two 
years. Thus, for a fixed computational time, if the computer 
cannot solve an instance when using a polynomial-time al- 
gorithm, one need not wait long as the exponential growth 
of computing power will reasonably quickly overtake the 
polynomially large cost. However, if the algorithm runs in 
exponential-time, one may be forced to wait several lifetimes 
in order for an instance to become soluble even if Moore's 
law continues indefinitely. Complicating matters, the expo- 
nential growth according to Moore's law is expected to cease 
sometime this century; hence the recent emphasis on quantum 
computation. 

The time complexity can be characterized using any com- 
putationally equivalent model. In the context of computer sci- 
ence, equivalent means that one model can simulate the other 
with an overhead that scales polynomially as a function of 
system size. These mappings respects the boundary between 
efficient and inefficient algorithms described earlier. Turing 
machines and circuits are two typical models discussed in the 
context of computational complexity. 

The Turing computer or Turing machine was introduced by 
Alan Turing 32 in 1936 before transistors and electronic cir- 
cuits and formalizes the idea of a computer as a person to 
whom computational instruction could be given. Turing! 3 -^ in- 
troduced the concept in order to answer David Hilbert's chal- 
lenge to decide if a polynomial has roots which are integers 
using a "finite number of operations." Turing proved that this 
was not possible using his Turing machine. An illustration de- 
picting the salient features of the Turing machine is given in 
Fig.0 

The circuit model of computation has been more widely 
used when generalizing to the quantum setting although the 
quantum Turing machine formulation does exist) 33 ! 34 ! The 
time complexity in the circuit model is characterized the num- 
ber of circuit elements (or gates) used to compute the ac- 
cept/reject boolean value corresponding to an input. So long 
as a universal set of gates, e.g. NAND, FAN-OUT, and FAN- 
IN for classical circuits, are used, the time complexity will 
only differ by a polynomial factor. The key caveat to consider 
when using this model is the notion of uniformity. To properly 
define an algorithm in the circuit model, a way of generating 
circuits for all input sizes must be given concisely. A circuit 
family for an algorithm is uniform if there is a polynomial- 
time algorithm that specifies the circuit given the input size. 

2.2 Computational problems and problem instances 

In computer science, computational problem often refers to 
decision problems which are collections of yes/no questions 
that an algorithm can decide. Each question is called a prob- 
lem instance and the answer is called the solution to that in- 
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Fig. 1 Illustration of a Turing machine that accepts input strings be- 
ginning with '01' or '10.' The computer is specified by the set of 
states (modes of operation), an alphabet of symbols for the scratch 
space, and the transition rules. The computer performing an al- 
gorithm reads the symbol on the current line, then based on the 
computer's current state and the given transition rules, the computer 
changes the symbol of the current line, changes its current state and 
moves either up or down. The computer halts when it enters the 
ACCEPT or REJECT states which indicate the output of the compu- 
tation. In the figure, the possible states are {BGN, CHK, SWP, RE- 
JECT, ACCEPT}, the alphabet is {'0','1'}, and the transition rules 
are listed in the box to the left. Depicted is the second step of a com- 
putation that began with input '10'. 



stance. Although we focus on decision problem, in many sit- 
uations the decision problem can be used to extract numerical 
values by asking sufficiently many yes/no questions. 

A collection of decision questions with affirmative answers 
such as, "Is 149 prime?" or "Is 79 prime?," is called a lan- 
guage. Languages define decision problems where only the 
accept instances are included in the language. If an algorithm 
accepts (returns "True" or outputs 1) on all strings contained 
in a language then the device is said to recognize the language. 
If the computer is also required to halt on all inputs, then it is 
said to decide the language. The difference being, a Turing 
machine may recognize a language but fail to halt on some 
instances, i.e. run forever. The problem of deciding if a Tur- 
ing machine algorithm halts is called the HALTING problem. 
This problem has a rich and storied history beginning with the 
first paper of Turing ^ . In Fig. [I] the algorithm decides the 
language of strings beginning with '01' or '10.' 

As an example from chemistry, the computational problem 
of deciding if a molecule has a dipole moment, DlPOLE, is 
a collection of questions: "Does x have a dipole moment?" 
Here, x is a string representing the molecule in each instance. 
The questions 'Does BCI3 have a dipole moment?', and 'Does 
NH3 have a dipole moment?' are instances of this computa- 
tional problem. The string x="BCi3" is not in the language 



Dipole since "BCI3" does not have dipole moment while 
"NH3" is in the language since it has a dipole moment. A rea- 
sonable modification of the problem would be 'Does x have a 
dipole moment greater than 0.1 Debyes?' so that small dipole 
moments may be ignored. 

It may also be necessary to promise that x represents a 
molecule in a specific way such that ill formatted inputs can be 
ignored. This is accomplished using promise problems where 
the promise would be that x is indeed a string that properly 
encodes a molecule. Promise problems also arise in order to 
account for issues of numerical or experimental precision as 
illustrated in the next paragraph. 

Promise problems also play a key role throughout the re- 
mainder of this text since we are discussing physical proper- 
ties where infinite precision is neither required nor expected. 
As an illustration, we use an example from thermochemistry, 
where language A is the set of strings corresponding to (ideal) 
gases with heat capacity at constant pressure, C„, less than 
some critical value. Imagine you have an unreasonable lab 
instructor, who gives you a substance that has heat capacity 
extremely close to the critical value that decides membership 
in A. It may take a large number of repetitions of the experi- 
mental protocol to be confident that the substance belongs or 
does not belong to language A. A reasonable lab instructor 
would announce at the beginning of lab that all the substances 
he is handing out for the experiment are promised to be at 
least one Joule per Kelvin away from the critical value. Given 
this promise, the student will be able to complete their lab in 
a single lab period. Without such a promise, it may take the 
student several lab periods to repeat the experiment in order to 
establish a sufficiently precise value of the heat capacity to de- 
cide the instance. Instead when using a promise, if the student 
has a compound that violates the promise then the instructor 
would give full credit for either answer. 

More formally, in the usual decision problems, the prob- 
lem is defined using one language, L, specifying the accept 
instances. The language of reject instances is all strings not in 
L. Promise problems differ in that they are specified using two 
languages: one for the accept instances, L accep[ , and a separate 
language for the reject instances, L re j ect . If an instance is not 
in either of these languages, then the promise is violated and 
the computation can terminate with any result. 

Reusing the DlPOLE example, if we can reformulate the 
problem to find out about non-zero dipole moments as a 
promise problem to account for experimental error. Now the 
computational task of DlPOLE is: 'Given a molecule, x, a 
trial value for the dipole dj, and an error tolerance 8, de- 
cide if the dipole moment of x is greater than dj + 8 or less 
than dj — 8, promised that the dipole moment is not between 
dj ± 8! In this case, there are two languages defining the 
problem L^^+g and L ( j <l i T _$. If a molecule has a dipole 
moment of exactly dj it would violate the promise and the 



Computer Instructions 

Start at line 1 in state BGN 
Transition Rules 

If in state BGN and in the current line there is a... 
'0', then change to state 'CHK', move down 
'1', then erase '1', write '0', change to state 'SWP', move down 

If in state SWP and in the current line there is a... 
'0'...then erase '0', write '1', change to state 'BGN', move up 
T...then erase '1', write '0', change to state 'BGN', move up 

If in state CHK and in the current line there is a... 
'0'... change to state 'REJECT', halt 
'1'... change to state 'ACCEPT, halt 

If line is blank ... change to state 'REJECT, halt 
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experimenter or computer does not need to answer or can re- 
spond with any answer. The value 8 allows us to meaningfully 
define problems in the presence of errors resulting from imper- 
fect experimental measurements or from the finite precision of 
a computing device. 

2.3 Computation reductions 

The idea of computational reduction is at the heart of classify- 
ing computational problems. Reducibility is a way of formal- 
izing the relationship between problems. Essentially, it asks, 
"If I can solve one problem, what other problems can I solve 
using the same resources?" 

There are two main types of reductions: Turing reductions 
and Karp reductions. For the reduction to be useful, they must 
be limited to polynomial time. The polynomial-time Turing 
reduction, also called the Cook reduction, of problem A to 
problem B uses solutions to multiple instances of B to decide 
the solutions for instances of A. The solutions to instances 
of B are provided by an oracle and each questions is called a 
query. Algorithms deciding A which require only polynomial 
queries of the oracle for B are efficient whenever the oracle for 
B answers in polynomial time. 

The other type of reduction is called a Karp reduction or 
polynomial transformation. If, instead of an oracle, there is an 
embedding of instances of problem A into instances of prob- 
lem B, then an efficient solution for all instances of problem 
B imply efficient solutions for instances of problem A. This 
transformation is a Karp reduction. When A can be reduced to 
B, under either Karp of Turing reductions, it is denoted A < B. 

To illustrate the difference, we use examples from thermo- 
dynamics. Consider trying to determine if the heat capacity at 
constant pressure, C„, of an ideal substance is less than some 
critical value (language A). Now, suppose an oracle, probably 
in the form of a bomb calorimeter, has been provided for de- 
ciding if the heat capacity at constant volume, Cy, of a given 
substance is above or below a critical value. Call language B 
the set of strings labeling substances with Cy below this value. 
By adding the number of moles times the ideal gas constant to 
the critical value, one can determine membership in language 
A via the formula C p = Cv + nR. Since each instance of A 
can be embedded into an instance of B, language A is Karp 
reducible to B and we can write A < B. 

Suppose instead, an oracle for evaluating if, at fixed pres- 
sure, the internal energy, U, of the given substance at a given 
temperature is less than some critical value (language C) is 
given. Then by evaluating the internal energy at two differ- 
ent temperatures, the heat capacity at constant pressure can be 
bounded by numerically estimating the derivative. Because 
the reduction has to use the oracle for C more than once to 
decide instances of A, language A is Turing reducible to C and 
A<C. 



2.4 Basic complexity classes 

Equipped with the basic definitions from computer science, 
we now introduce the concept of computational complex- 
ity classes. This will give some insights into why quantum 
chemistry is difficult. We will introduce six basic complex- 
ity classes to classify the time complexity of decision prob- 
lems. The first three complexity classes are characterized by 
algorithms that can decide instances in time proportional to a 
polynomial of the input size; the other classes are character- 
ized by polynomial time verification of problem instances. 

In table[T] the three polynomial time complexity classes are 
listed. First, P, is the class of all decision problems with in- 
stances that can be accepted or rejected by a Turing machine 
in polynomial time. If the Turing machine has access to an 
unbiased random number generator, then the complexity class 
of problems that can be decided in polynomial time is called 
BPP. The term "bounded error" refers to the requirement that 
the probability of acceptance and of rejection must be bounded 
away from half so that repetition can be employed to boost 
the confidence in the answer. An important class of problems 
falling into this complexity class are efficient Monte Carlo 
simulations often grouped under the umbrella term quantum 
Monte CarlcJ^JtlH use( j f or com p U ting electronic structure in 
chemistry. Lastly, the complexity class BQP is character- 
ized by problems soluble in polynomial time with quantum 
resources. In the quantum computational models, the quan- 
tum algorithm is conceptually simpler to think of as a unitary 
circuit, U, composed of unitary circuit element that affect, at 
most, two quantum bits. As mentioned earlier, the number of 
gates used determines the time complexity. The outcome of 
the algorithm with input \Input) is "accept" with probability 
| {Accept \U\Input) \ 2 . Similarly, for the "reject" cases. 

So far, the discussion has centered on complexity classes 
containing problems that are solvable in polynomial time 
given access to various resources: a computer (P), a random 
number generator (BPP), and a controllable quantum system 
(BQP). Now, our attention turns to the class of problems that 
can be computed non-deterministically. The original notion 
of non-determinism is a Turing machine whose transition rule 
maps the computer state and the tape symbol to any number of 
possible outputs. In essence, this is a computer that can clone 
itself at will to pursue all options at once. NP is the class of 
problems that could be solved in polynomial time by such a 
computer. Whether a deterministic polynomial time computer 
can be used to simulate such a computer is a restatement of the 
famous open question in computer science: "Does P = NP?" 
This question was selected as a millennium problem by the 
Clay Mathematics Institute which has offered a one million 
dollar prize for a correct proof of the answer. 

Rather than resorting to imaginary computers, the non- 
deterministic classes can be defined by considering verifi- 
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Fig. 2 An example of the non-deterministic complexity class QMA. Non-deterministic problems can be thought of as games. In the Merlin- 
Arthur games, the proof verifying the validity of input x is magically (hence non-deterministically) given by a wizard Merlin who is prone to 
deception. The verifier, called Arthur, is trying to decide if he should accept or reject input x. When Merlin gives Arthur a valid proof that 
verifies x should be accepted, Arthur should accept with high probability (completeness) and when Merlin gives him an proof of something 
incorrect, he should be able to spot the inaccuracy (soundness). In the drawing since Arthur has access to quantum bits (qubits), the complexity 
class depicted is QMA. Figure composed from images found on Wikipedia under the Creative Commons License. 



cation of solutions which have been obtained in some non- 
deterministic way, see Fig. [2] For every "accept" instance, x, 
there exists at least one proof state y such that the verifier re- 
turns "accept" in polynomial time. If x is should be rejected, 
then for every proof state y, the verifier should output "re- 
ject", again, in a polynomial amount of time. Each of the 
non-deterministic classes listed in table [T] are characterized 
by the computational power of the verifier. Note that P is a 
subset of NP because any problem that can be easily solved 
can be easily verified. 

2.5 Completeness and hardness 

The goal of classifying computational problems into complex- 
ity classes motivates introduction of the terms hard and com- 
plete. The classification of a computational problems as hard 
for a complexity class means if an algorithm can solve this 
problem, then all problems in the class can be solved effi- 
ciently. In other words, this problem contains all problems in 
the class either via Karp or Turing reductions. More precisely, 
a language B is hard for a complexity class CC if every lan- 
guage, A, in CC is polynomial-time reducible to B, i.e. A < B. 

Complete problems for a class are, in some sense, the hard- 
est problems in that class. These problems are both a member 
of the class and are also hard for the class. That is, a language 
B is complete for a complexity class CC if B is in the complex- 
ity class CC and for every other language A in CC, A < B. 

The simplest illustration of the difference between hard 
and complete computational problems is the distinction be- 
tween optimization problems and decision problems. Opti- 
mization problems such as finding saddle points or minima 
in energy landscapes are frequently encountered in computa- 
tional chemistry, however these problems are not in the com- 
plexity class NP. However, if one can perform optimization 



quickly, then responding with the answer to yes/no questions 
about the solution would be simple. Thus, optimization prob- 
lems can be classified as NP-hard but not NP-complete since 
the computational task of optimization is, in a sense, harder 
than just answering yes or no. 

3 Electronic structure and other Hamiltonian 
problems 

Now, armed with the key ideas from computer science, we 
return to our original inquiry into why quantum chemistry 
is hard. We answer this questions by exploring the compu- 
tational complexity of three different widely used methods 
for computing electronic energy in computational chemistry: 
Hartree-Fock, two-electron reduced density matrix methods, 
and density functional theory. 

Before delving into the specific complexities of these prob- 
lems, we first establish notation by reviewing necessary con- 
cepts from quantum chemistry, then discussing computational 
problems concerning classical and quantum spin Hamiltoni- 
ans. We end the section with mappings between systems of 
electrons and spin systems. 

3.1 Quantum chemistry and second quantization 

In quantum chemistry, the annihilation, {a^}, and creation op- 
erators, correspond respectively to removing and adding 
an electron into one of M single particle wave functions, 
{^(xi)}^. The single-particle wave functions are called or- 
bitals and the set of orbitals is typically called the basis set. 
To include the electron spin, the single particle orbitals 0^ are 
functions of spatial coordinates and a spin variable which are 
collectively denoted x. 
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Table 1 Polynomial time complexity classes are separated by the re- 
sources the computer has access to while the non-deterministic poly- 
nomial time complexity classes are characterize by the resources of 
polynomial time computational verifiers 



Class 


Name 


Computer type 
taking only poly, time 


1 1 

r 


Polynomial Time 


Turing machine 


BPP 


Bounded Error 
Probabilistic Polynomial 
Time 


Turing machine with access 
to random number generator 




Bounded Error 
Quantum Polynomial 
Time 


Turing machine with access 
to quantum resources 


Class 


Name 


Verifier's computer type 
taking only poly, time 


NP 


Non-deterministic 
Polynomial Time 


Turing machine 


MA 


Merlin-Arthur 


Turing machine with access 
to random number generator 


QMA 


Quantum 


Turing machine with access 


Merlin-Arthur 


to quantum resources 



The electronic spin will play an important role when dis- 
cussing the connections between systems of electrons and sys- 
tems of quantum spins. Since electrons are spin-j particles, 
the electron spin is either up or down. Accordingly, there are 
M/2 orbitals with spin up and M/2 with spin down. Unless 
explicitly noted, the summation over orbitals includes a sum- 
mation over the spatial and spin indices of the orbitals. 

Anti-symmetry of the Af -electron is enforced by the canon- 
ical anti-commutation relations for fermions, i.e. electrons, 



l pj a q\ + 



= 0, 



l Pi a q\ + 



= 8 pq l. 



(1) 



In chemistry, typically the electronic structure is the pri- 
mary concern, but if interested in vibrational structure, we 
would have to consider bosonic canonical commutation re- 
lations: [b p ,b q ] = b p b q ~b q b p = and [b p ,b q ] = S pq l. The 
vacuum state, \vac), is the normalized state that the cre- 
ation/annihilation operators act on and it represents the system 
with no particles. 

Acting on the vacuum state with a strings of N distinct cre- 
ation operators yields A^-electron wave functions. The most 
general A^-electron state within a basis set is, 



hi 



I*) = H CKa k l a K 2 --- a K N \ Vtl > 
K 



(2) 



with K = (K\ ,Kn) and the complex valued Cg are con- 



strained such that La^Cr - ! 2 = 1. To convert between cre- 
ation and annihilation operators and the coordinate repre- 
sentation, one uses field operators: 0(x) = Y.k^k{^) a k as m 
*P(xi,--- ,x N ) = (vac\^(x l )---^{x^)\'¥}/VNl. The factor of 
AH accounts for permutation symmetries in the summation. 
Each state or I 1 !') is called an A^-electron pure state. A valid 
Af-electron mixed state, 



,(A0 



(3) 



has Y,Pi = 1 where each j*?;) is an A^-electron wave function. 
The k-th order reduced density matrix, abbreviated &-RDM, is 
defined using the field operators: 



1 



Tr 



(4) 
(5) 



With (N)^ = N\/ (N — k)\, the reduced density matrices are 
normalized to unity and, in equation |5]), the trace sums over 
the expectation values of states from a complete set of N — k- 
electron wave functions. 

Lastly, we define the primary computational chemistry 
problem encountered in computation of electronic structure. 

Computational Problem: ELECTRONIC STRUCTURE. The 
inputs are the number of electrons, N, a set of M orbitals, a 
static configuration of nuclei, a trial energy, Ej, an error toler- 
ance 8 and Hamiltonian, 



H, 



elec 



eN : 



Vf, N )aUi 



T tPP t t 

■V ijkl a\a)a k ai 
(6) 

with T e the electronic kinetic energy operator, V ee the electron- 
electron interaction operator and V e N the electron-nuclear in- 
teraction. The task is to decide if the ground state energy is 
less than Et — S or greater than Et + 8 promised that the en- 
ergy is not between £Y ± <5. 

Methods related to getting approximate solutions to ELEC- 
TRONIC Structure have complexities ranging from NP- 
complete to QMA-complete as shown in sections [4] and [5] 
The complexity of Electronic Structure itself is not yet 
clear. In the literature, only Hamiltonians with more flexibil- 
ity, namely a local magnetic field, have been demonstrated to 



be QMA-complete as shown see section 6.1 



3.2 Classical and quantum spin Hamiltonians 

The problem of estimating the ground-state energy of Hamil- 
tonians of different forms lies at the intersection of computer 
science, physics and chemistry. In this section, we define 
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computational problems related to both classical and quantum 
spin Hamiltonians. Spin systems have lon g bee n known to 
provide fertile ground for complexity theoryG^SEl Given a 
Hamiltonian, deciding if a there exists a spin configuration, 
e.g. (ftl • • • )> tnat satisfies a certain property, e.g. has energy 
below a certain threshold, can be difficult even when check- 
ing the property for each configuration is easy. Computations 
designed to decide properties of this sort are often in the com- 
plexity class NP. If, additionally, that particular property can 
embed instances of any other NP properties, then that property 
corresponds to an NP-complete computational problem. 

Computational Problem: IsiNG. The inputs are N classi- 
cal spins with possible values ±1, trial energy, Ej, and Ising 
Hamiltonian 

neighbors 

Hising = — ^ J'jSiSj (7) 

(y) 

having 7, ; either or ±1. The task is to decide if there is a 
configuration of spins c = (s\,S2,- ■ • ,sn) such that the energy 
of the Ising Hamiltonian is less that Ej or if all configurations 
have energy above Ej. 

Unlike the other Hamiltonian problems discussed in this re- 
view, IsiNG does not require a promise on the precision be- 
cause the energy values of the Hamiltonian are integer-valued 
instead of continuous. 

Although some instances of the Ising lattice can be solved 
analytically, e.g. one-dimensional Ising lattices or the solu- 
tion by Onsager^lfor two-dimensional Ising lattices with uni- 
form couplings, Barahona 39 showed that problem ISING is 
NP-hard when interactions are restricted to nearest neighbors 
on a 2 x L x L lattice. Since this model is only of tangential 
interest to electronic structure, details of this proof are omit- 
ted. Since the "accept" instances can be verified in polyno- 
mial time given a configuration c with energy less than Ej, 
this problem is in the complexity class NP and hence, is NP- 
complete. Further results show that all non-planar lattices, 
even in two-dimensions, are NP-completeP. 

To provide examples of Hamiltonian problems that are 
complete for the quantum analogue of NP, we turn to the quan- 
tum analogues of the Ising Hamiltonian: quantum spin Hamil- 
tonians. To define similar quantum Hamiltonians with QMA- 
complete properties, it turns out that one only needs spins with 
angular momentum of i which are often called qubits. The in- 
teractions are now expressed using tensor products of the Pauli 
sigma matrices with matrix representations 



" 1 " 




" -i ' 


, o 7 - = 


' 1 " 


1 




i 


-1 



Additionally, a is the identity matrix. 



A widely studied class of qubit Hamiltonian problems 
whose computational complexity was first investigated by Ki- 
taev et al.^ is the problem of computing, to polynomial ac- 
curacy, the ground state energy of fc-spin local Hamiltonians. 
This set of problems is referred to as A:-Local HAMILTO- 
NIAN in the computational complexity literature followingED 
but we use the name £-Local-Spin Hamiltonian to em- 
phasize the nature of the Hamiltonian. 

Computational Problem: &-L.OCAL-SPIN HAMILTONIAN. 
Given A:-spin-local Hamiltonian acting on N spins, 

m 

H kQM A = - E Jcpotl®o%® (9) 

D=(d h ...,d k ) 
C=[c\ c k ) 

where dj £ {x,y,z,Q}, c, labels a spin, there are m =poly(A/) 
terms, and \ Jc,d\ < 1, decide if the ground state energy is less 
that Eo — 8 or if the ground state energy is greater than Eq + 8 
with 8 < 1 /poly (N) promised that is is not between Eq ± 8. 

Note that the problem is defined so that the error scales rela- 
tive to the number of non-zero terms m in equation (|9j. To un- 
derstand this, consider the following set of instances with er- 
ror tolerance 8 = 0.5. Consider Hamiltonian, H, with ground 
state energy 17.7 and Hamiltonian, H' , with energy 17.6. As 
it stands, the two Hamiltonians cannot be distinguished. How- 
ever, by including each term in H and H' ten additional times, 
the ground state energies are now 177 and 176, which can be 
resolved at error tolerance 8. By considering m as a polyno- 
mial in N, one cannot indefinitely rescale the energy to effec- 
tively shrink the error tolerance. 

The first demonstration of a QMA-complete problem^] re- 
quired k = 5. Subsequently k was reduced^ to 2. Finally, 
the problem was shown^ to remain QMA-complete even 
when Jc.d is non-zero only if two spins are spatially adja- 
cent on a 2D lattice. The problem can be further restricted 
by including only a limited set of two-spin interactions and 
still remain QMA-complete 45 . Complexity results concern- 
ing qubit Hamiltonians as well as other variants with higher 
dimensional spins and more restrictive lattices were recently 
reviewed by OsborneG2 

3.3 Relationships between spin systems and electronic 
systems 

Since chemists are not necessarily interested in qubit or quan- 
tum spin systems, we now discuss connections to systems in- 
distinguishable electrons. In this subsection, we will examine 
how to embed a system of spins into an electronic Hamilto- 
nian and how to embed electronic systems into spin Hamilto- 
nians. Lastly, we use these connections to give an example of 
a fermionic system (although not ELECTRONIC STRUCTURE) 
that is QMA-complete. 
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Given an electronic Hamiltonian, there is a orbital pair 
pseudo-spin mapping^ES that is used to embed qubit mod- 
els to the ground state of the electronic systems. To em- 
bed N spins, there must be M = 2N orbitals and N electrons, 
i.e. half-filling, and each quantum spin is identified with a pair 
of fermionic modes. The embedding requires translating each 
spin \qi) = a\ T,-) +]8| |,) to fermionic operators: a aj^ + j3 atj_. 
The Pauli operators appearing in the spin Hamiltonian, then 
become single fermion terms, e.g. h[ja\aj. As important ex- 
amples consider, 

of =Ui)<t;| + lti)Ui| 4 fl ''t + fl W (10a) 

of = i(U»Xti|-|ti>Ui|) /(«V, t -4«4) (10b) 

of =|ti> (til- \U)(U I fl^t-ajfljj. (10c) 

Taking tensor products of the Pauli matrices leads to two- 
fermion terms, e.g. hij^a^a^a^ai. The final concern is pre- 
venting double occupancy within a pair of sites which would 
invalidate the pseudo-spin interpretation. This is 

handled™] 

by including an additional term which penalizes invalid elec- 
tronic configurations. For each pair of electronic modes, the 
two-fermion penalty 

Pi = C a'^a^a]^ (11) 

is added to the electronic Hamiltonian. Because penalty P = 
Y,i Pi commutes with the fermionic Pauli matrices, the ground 
state still corresponds to the solution of the spin Hamiltonian. 
The constant C can be selected as a low order polynomial in 
the system size to ensure that the ground state remains in a 
valid pseudo-spin stateQ Note that we have not explicitly re- 
lied on the anti-commutation properties and a nearly identical 
construction exists for bosonic systems^. 

A second connection is given using well established tech- 
niques developed to translate certain spin systems to sim- 
pler non-interacting fermionic systems that can be exactly 
The Jordan- Wigner transform 51 provides this 
connection by mapping fermions to spin operators such that 
the canonical anti-commutation relations are preserved. The 
Jordan-Wigner transform between N fermionic creation and 
annihilation operators and the Pauli matrices acting on N spins 
is given by 

aj & l^- 1 (7+ <g> a z ® N ~ 3 (12a) 
at^l^-^a-fSxr^-j (12b) 

where o+ = = UXtl and a = = | t>U I- The 

qubit state | t • • • t) corresponds to the vacuum state and the 

tMore precisely, the norm of the Hamiltonian \\H\\ is upper bounded by the 
sum of the individual terms. Since there are, at most, O(N) 4 terms the norm 
of the total Hamiltonian must scale less than a fourth order polynomial in the 
system size. 



string of a z operators, preserve the commutation relations in 
equation ([I} since <J Z and ff* anti-commute. More sophisti- 
cated generalizations of the Jordan-Wigner transform, reduce 
the number and the spatial extent of the spin-spin interac- 
tions! 54 " 56 ! 

The orbital pair pseudo-spin mapping and the Jordan- 
Wigner transformation complement each other and will be 
used repeatedly throughout the remainder of the article. The 
orbital pair pseudo-spin mapping requires carefully engi- 
neered penalties and fixes the total number of orbitals. Thus, it 
is primarily useful when translating spin systems to electronic 
systems. By contrast, the Jordan-Wigner transformation is pri- 
marily useful in the other direction, that is, when translating an 
arbitrary electronic system to a quantum spin system. The next 
subsection gives a concrete example of how these connections 
are employed when studying the computational complexity of 
electronic systems. 

3.4 Generic local fermionic problems are QMA-complete 

The orbital pair construction allows one to immediately ascer- 
tain that the ground state energy decision problem for Hamil- 
tonians containing two-fermion interactions, 

H 2f = Y,kj<4 a j+hjklOi a ) a k a h (13) 

is QMA-hard 48 : This follows via a Karp reduction using 
the orbital pair pseudo-spins as in equation ( fT0] >. Since only 
two-spin interactions in equation |9]) are required for QMA- 
completeness 43 , each spin-spin interaction term translates to 
two-fermion terms under the pseudo-spin mapping. Since the 
pseudo-spin construction can be extended to bosonic systems, 
the two-boson ground state energy decision problem is also 
QMA-harcPl 

The ground state energy of the electronic state can be ver- 
ified to be above or below Ej ± 8 by a BQP quantum com- 
puter given the proof state. This implies that the problem is 
QMA-complete. We sketch the idea relying on well kno wn re- 
sults about BQP quantum simulation of chemical systems^EIl 
The Jordan-Wigner transformation is used to translate the two- 
electron Hamiltonian into a qubit Hamiltonian. If Merlin pro- 
vides the fermionic ground stat^*] the energy of correspond- 
ing qubit state can be determined through simulation ! 8 ! 57 ! 58 ! 
of the evolution under the qubit version of the two-electron 
Hamiltonian. This quantum evolution can be simulated effi- 
ciently on a quantum computer 59 ! 60 ! The resulting evolution 
in the time-domain is Fourier transformed to extract the en- 
ergy of the ground stat d 5 ! 8 ! 11 ^, allowing the verifier to deter- 
mine whether to accept or reject the instance. This argument 
carries through for the bosonic case as well using the bosonic 
equivalent of the Jordan-Wigner transform^Sl. 



| The verifier requires multiple copies of the state to ensure that the state has 
exactly N electrons. 
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4 Hartree-Fock 



4.1 Hartree-Fock is NP-complete 



Hartree-Fock is one of the most important quantum chemistry 
techniques as it typically recovers about 99% of the total elec- 
tronic energy. Hartree-Fock is known to be a weak approxi- 
mation in many instances, but it is the basis for more sophis- 
ticated (post-Hartree-Fock) methods which improve upon the 
Hartree-Fock wave function. Furthermore, Hartree-Fock pro- 
vides the mathematical framework for the widely adopted no- 
tion of molecular orbitals used throughout chemistry. 

The implementation of the Hartree-Fock algorithm requires 
evaluating and manipulating 0(M 4 ) two-electron integrals. 
When this computation dominates the runtime, the algorithm 
scales as somewhere between 0(M 2 ) and (9(M 3 ). However, 

since Hartree-Fock solves a nonlinear eigenvalue equation 

m 

through an iterative method , the convergence of self con- 
sistent implementations is the key obstacle that prevents the 
worst case scaling from being polynomial. The computational 
complexity result proving that the worst case scaling cannot 
be polynomial unless P=NP was provided in an unpublished 
appendix of Schuch and Verstraete 49 available on the arXiv 
preprint server. 

For the purposes of this article, the Hartree-Fock procedure 
can be succinctly explained as the minimization of the energy 
of an A^-electron system given M basis functions with the re- 
striction that in the expansion found in equation |2]), all Ck are 
zero except one. Explicitly, 



E„ F = min (^\H elec \W) 

*=n,(^)*'|v«c>, 



(14) 



Here, the single Fock state corresponding to the minimal value 
of Ehf is called the Hartree-Fock state: ^hf- The optimized 
set of creation operators, 



(15) 



place and remove electrons from the set of molecular or- 
bitals, V}( x ) — H,f Cij<f>i(x). The formal computational prob- 
lem Hartree-Fock will be defined analogous to the other 
Hamiltonian problems with a promise given to account for 
precision. 



To prove that any other NP problem can be mapped to the 
HF problem (NP-hardness), we can use a Kara reduction to 
embed instances of IsiNG into instances of HP 4 -^. Using the 



fermionic version of a z given in equation ( 10 1, Hi s i ng only re- 
quires two-electron interactions, 



SiSj 



(16) 



The satisfying assignment of the IsiNG instance is some spin 
configuration \s\S2 ■ • • %) f° r N = 21? with Sj as either j" or \.. 
The correct and exact Hartree-Fock solution for the iV-electron 



wave function should assign bj — at when Sj =t and b] — 

ajj_ when s, =|. Thus, Ising<Hartree-Fock under a Karp 
reduction. 

The promise on the error tolerance, 8, given in the specifi- 
cation of the problem is necessary because of the pseudo-spin 



representation used in equation (16i. Corrections to the en- 
ergy due to the pseudo-spin representation can be computed 
using perturbation theory with the unperturbed Hamiltonian 
given by equation (jTTJ and the converted Ising Hamiltonian, 
c.f. equation ( 16 1, as the perturbation. The first order correc- 
tions in small parameter C~ l are the Ising energies and the 
errors due to the pseudo-spin mapping arise at second order in 
C _1 . A coarse estimate for C is obtained by multiplying the 
number of non-zero terms by the maximum absolute value of a 



coefficient occurring in equation ( 13 1. Since we map the prob- 
lem to ISING, there are 0(N 2 ) terms in the summation and the 
maximum absolute value of each term is unity. Hence, C is 
estimated as 0(N 2 ). So long as 8 < 0(N~ 2 ), the first order 
corrections occurring at order C _1 can be distinguished and 
the ground state Ising energies can be recovered. 

To prove inclusion of HF in the complexity class NP, an 
algorithm for verifying the energy in polynomial time must 
be given. If the coefficient matrix C describing the correct 
orbital rotation is given, then the energy is calculated easily 
using Slater-Condon rules 62 . Thus, given the answer, the ac- 
cept/reject conditions are quickly verified. Since the prob- 
lem can also be quickly verified, HARTREE-FOCK is NP- 
complete. 



Computational Problem: HARTREE-FOCK. The inputs are 
the number of electrons, N, a set of M orbitals, a trial energy, 
Et, an error tolerance, 8 < 1 /poly (AT), and a two-electron 
Hamiltonian, c.f. equation ( fT3| . The sum of the absolute val- 
ues of the coefficients is required to scale less than a poly- 
nomial in N. The task is to decide if the Hartree-Fock en- 
ergy Ehf = (^hf {H^hf) is less than Et — 8 or greater than 
Et + 8 promised that E^f is not between Et ± 8. 



5 2-RDM methods 

Many computational chemistry algorithms seek to minimize 
the energy or other properties by manipulating the wave func- 
tion but to quote Coulson^l, "wave functions tell us more than 
we need to know. . . All the necessary information required for 
energy and calculating properties of molecules is embodied in 
the first and second order density matrices." Extensive work 
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has been done to transform this remar k into a host of computa- 
tional methods in quantum chemistry! 64 | 65 | j n t jy S sec ti on anc j 
the following, we review the prominent computational com- 
plexity results related to this body of work. 

With respect to the 2-RDMs, it is easy to evaluate properties 
of the electronic system since the Hamiltonian only contains 
two-electron terms. However, difficulties arise when deter- 
mining if the 2-RDM is valid or invalid. Unconstrained opti- 
mization of the energy using 2-RDM can lead to ground state 
energies of negative infinity if the validity of each 2-RDM can- 
not be determined. While the criteria for validity have recently 
been developed^!, deciding the validity of each 2-RDM us- 
ing such criteria is known as the N -representability problem. 
Mazziotfr^l also provides a proof that the following problem, 
jV-Rep, is at least NP-hard. In the next subsection, we fol- 
low Liu et al.^ to demonstrate the stronger conclusion that 
Af-REP is QMA-complete. 

Computational Problem: A^-Rep. The inputs are the num- 
ber of electrons, N, an error tolerance j3 > l/poly(A r ), and a 
2-RDM, n^ 2 \ The task is to decide (i) if /J.^ is consistent with 
some Af-electron mixed state, p( N \ or (ii) if jx^ is bounded 
away from all 2-RDMs that are consistent with an A^-electron 
state by at least j3 [^promised that either (i) or (ii) is true. 

5.1 N-Rep is QMA-complete 

To show QMA-hardness, a Turing reduction is used to show 2- 
Local-Spin Hamiltonian< N-Rep. Before proceeding to 
the Turing reduction, note, since A^-Rep only considers a fixed 
number of electrons, the one-electron operators are not needed 

a&a\aj=a] (xf a\aA aj/(N-l). Each valid 2-RDM can be 

represented by an (¥) dimensional matrix with (^f) 2 — 1 
independent parameters'] A complete set of observables, such 
as the projection onto each matrix elemenjj] is then used to 
characterize the space of 2-RDMs. 

The space of valid 2-RDMs is convex; that is if 

(2) (2) (2) 

jUj , H2 , • • • , Ml are valid 2-RDMs, then the convex sum 

Zjij^ with Ljnj = 1 is also a valid 2-RDM. This follows 
directly from the convexity properties of sets of density matri- 
ces and probability distributions. 

With access to an oracle for A^-Rep, the boundaries of the 
convex set of valid 2-RDMs can be characterized. Because the 



§The appropriate metric is the trace distance, d,,-(A,B) = \\A — B\\ tr = 
jTr[y/ (A — BY(A — B)], which generalizes the distance metric from standard 
probability theory. 

1 Since each 2-RDM is hermitian, the complex off-diagonal elements occur in 
pairs and the diagonal elements must be real. Since the trace is normalized to 

unity, there a reduction of one degree of freedom, leaving (^) — 1 indepen- 
dent parameters. 

|| Liu et a/.ESl used a different set observables inspired by the Pauli matrices 
with more convenient mathematical properties. 



oracle is used multiple times throughout the verification pro- 
cedure, this is a Turing reduction. Since convex minimization 
problems can be solved efficiently and reliably, see e.g.^l the 
QMA-hardness is nearly demonstrated. The last point of con- 
cern addressed by Liu et a/. EH, are the errors introduced by 
the promise given on the oracle. To demonstrate that the algo- 
rithm remains robust in spite of such errors, the authors use a 
tailored version of the shallow-cut ellipsoid convex optimiza- 
tion technique^. 

Let us remark that the required error tolerance follows since 
the reduction relies on the QMA-completeness of 2-LOCAL- 
Spin Hamiltonian which includes a promise on allowed the 
error tolerance. 

To demonstrate that A^-Rep is QMA-complete, what re- 
mains is demonstrating that A^-Rep is in the complexity class 
QMA. If /i' 2 ' is a valid 2-RDM, then, relying on the Jordan- 
Wigner transform introduced earlier, Merlin can send polyno- 
mial copies of the correct A^-electron state, The veri- 
fier, Arthur, first checks that the number of electrons in the 
given state is N by evaluating Tr[/iW£ a\a^\. Then Arthur 
randomly picks observables from the complete set described 
earlier and tests that the expectation value of the state juW 
and jii' 2 ' match until convinced. Merlin might try to cheat 
by sending entangled copies of fjW but it can be proven that 
Arthur cannot be fooled based on the Markov inequalit)|**] 
See Aharonov and Regev^ for a proof of this fact. 

5.2 Restriction to pure states 

The pure state restriction of A^-Rep has interesting conse- 
quences. Consider, Pure A^-Rep, where the question is now: 
"Did the 2-RDM come from a A^-electron pure state (up to er- 
ror /3)?" This contrasts with the original problem where the 
consistency questions refers to an A^-electron mixed state. 

This problem is no longer in QMA since the verifier, Arthur, 
cannot easily check that the state is pure if the prover, Merlin, 
cheats by sending entangled copies of l 1 ?^). If instead, two 
independent unentangled provers send the proof state, the state 
from the first prover can be used to verify thepurity. The 
purity is checked using the pairwise swap testPH, a quantum 
algorithmic version of the Hong-Ou-Mandel effect in quantum 
optics^, on each of the supposedly unentangled states. When 
the given state is not a product state (or nearly so), Arthur will 
detect it. If the set of states from the first prover pass the test, 
the second set can be used to verify that randomly selected 
expectation values of |\pW) and )J,( 2 > match. The complexity 
class where two verifiers are used instead of one verifier is 
QMA(2). 



** Given any random variable X, the Markov inequality states (|X|) > aPr(|X| > 
a) with (X) defining the expectation value of random variable X and \X\ defin- 
ing the absolute value of X. 
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6 Density functional theory 

In this section, a further reduced description is examined: den- 
sity functional theory. Density functional theory (DFT) meth- 
ods are of profound importance in computational chemistry 
due to their speed and reasonable accuracy ^ 65 ! 73 !. In this sec- 
tion, the complexity of the difficult aspects of DFT is shown 
to be QMA-complete. This was first conjectured by Rassolov 
and Garashchuk 20 and rigorously proven by Schuch and Ver- 
straete^Sl. 

In DFT, the wave function is replaced by the one-electron 
probability density, n(x) = pW(x,x). The use of an object 
in three spatial dimensions to replace an object in 3N spatial 
dimensions without losing any information seems almost ab- 
surd, but the theoretical foundations of DFT a re well estab- 
lished™. The Hohenberg-Kohn fheoremEEHproves that the 
probability density obtained from the ground state wave func- 
tion of electrons in a molecular system is in one-to-one corre- 
spondence with the external potential usually arising from the 
static nuclear charges. Therefore, all properties of the system 
are determined by the one-electron probability density. 

In the proof of the Hohenberg-Kohn theorem, one of the 
most important, and, elusive, functionals of the one-electron 
density is encountered: the universal functional. In all elec- 
tronic Hamiltonians, the electrons possess kinetic energy and 
the electrons interact via the Coulomb interaction; only V e N 
due to the nuclear configuration changes from system to sys- 
tem. Separating these parts allows one to define the universal 
functional of DFT as 

F[n(x)]= min Tr[(T e +V ee )pW], (17) 

pW->n(s) 

which takes as input the probability density and returns the 
lowest possible energy of T e + V ee consistent with the prob- 
ability density. The nuclear potential energy can be directly 
determined efficiently using / dx n(x)V e N(x) using Gaussian 
quadratures or Monte Carlo sampling. Regardless of the 
method used to evaluate the integral, the domain does not (ex- 
plicitly) increase with the number of electrons. As a decision 
problem, we have the following: 

Computational Problem: UNIVERSAL FUNCTIONAL. The 
inputs are the number of electrons, N, a probability density 
«(x), a trial energy Ej and an error tolerance 8 < 1 /poly(7Y). 
The task is to decide if F[n(x)] is greater than Ej + 8 or less 
than Ej — 8 promised that F[n(x)] is not between Ej ± 8. It is 
required that the summation over the Hamiltonian coefficients 

i \ v i e fki\ + \ T i e j\ + \Vtf\ + \ V !T S \ scale less than P ol y(^)- 



2QMA 



heisenberg 



hubbard 



< 



6.1 Universal Functional is QMA-complete 

The demonstration that UNIVERSAL FUNCTIONAL is QMA- 
hard proceeds via a series of reductions. Specifically, Schuch 



and VerstraeteESI show that H 
H e iec where H\ < H2 means that instances of the ground state 
decision problem for Hamiltonian H\ can be embedded (Karp 
reduced) into ground state problem instances of Hamiltonian 
H2. In their proof, the authors utilize the magnetic field to en- 
code the problem instances causing H e / ec as defined in equa- 
tion (|6| to include an additional local magnetic field, V mag . 
Note that the local field only affects the spin of the elec- 
tron and does not require spin-dependent density functionals 
since the charge and spin do not couple. The Hamiltonian, 
Hiqma, was listed earlier in equation Again since the 
reduction relies on the QMA-completeness of the 2-LOCAL- 
Spin Hamiltonian problem, the promised error tolerance 
and the upper-bound of the coefficients of the Hamiltonian are 
required for QMA-completeness. 

The Hhubbard and Hheisenberg Hamiltonians are commonly 
encountered models in condensed matter physics and are of 
the form 



H, 



hubbard 



M/2 M/2 

E E ta h a u+ E Ua \\ a ]i a a a i\ 

M/2 

"EE Bfaf (18) 

i de{x,y,z} 

M M 



H, 



heisenberg 



E Z^M+L^f 



(19) 



de{x,y,z} ij 



The first Hamiltonian describes an electronic system where 
the Pauli matrices, a d , are expressed using orbital pairs as 
while the second Hamiltonian describes a 



10) 



in equation 

system of quantum spin. In both Hamiltonians, different prob- 
lem instances are embedded using the local magnetic field. 
The embedding of Hf, e i sen f, erg instances into Hhubbard follows 
directly from the orbital pair mapping described earlier in sec- 
The remaining two Karp reductions are more in- 



3.3 



tion 
volved. 

First, let us consider the reduction Hiqma < Hheisenberg- This 
Karp reduction follows along the same lines used in reduc- 
ing kQMA from k = 5 to k = 2 based on perturbative gad- 
gef^ 43 l 44 l widely used in quantum complexity proofs. A medi- 
ator spin splits the state space of the system into a low en- 
ergy and a high energy sectors. In the low energy sector, 
the perturbative coupling to the high energy states "mediates" 
new interactions in the low energy sector. For example, with 
H m = B m \(p m ) (<p,„| acting on the mediator spin and a perturba- 
tion V = <?f<J* <Jj + of CT^oJ, in the low energy space of spin 
m, there is, at second order, an effective interaction: of oj. 
Since the Hamiltonian problems refer to the ground state en- 
ergy, the high energy sector is not important. 

The reduction H 2 qma < H heisenberg requires: 

1. converting arbitrary strength couplings to constant 
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strength couplings: /i 2 of of H> /(of + o^of ); 
2. converting inequivalent couplings to equivalent cou- 



plings: of of h> of o^ + a^a 



B^B. 



3. and finally converting equivalent couplings to Heisenberg 
interactions: of o% H> £ d vftfl+Ld <^m a \- 

The full transformation requires 15 mediator spins where the 
local field on each of the mediator spins splits the system into 
low and high energy sectors. The specific local field depends 
on the interaction desired. The strength of the local field ap- 
plied to the mediator spin \B m | is usually very large to ensure 
that perturbation theory applies. So far this has limited the 
practical relevance of these constructions; for instance, field 
strengths required on the final set of mediator spins scales at 
nearly the 100th power of the system size. 

The remaining reduction, Hbubbard < H e [ ec , was heuristi- 
cally known since the Hubbard model phenomenologically 
describes electrons in solid state systems. Schuch and Ver- 
straete 49 rigorously demonstrate this reduction by careful ac- 
counting for error terms. They begin with a simple model used 
for studying solids where non-interacting (spin-less) electrons 
are subjected to a periodic delta function potential with M/2 
such sites. The orbitals of this system can be solved exactly. 
The resulting second quantized Hamiltonian has uniform ki- 
netic hopping terms: T^Uj) ta-aj. After incorporating the elec- 
tron spin, electron-electron interactions are introduced. The 
strength of this interaction is rescaled by changing the spatial 
distance between neighboring sites until only the electrons at 
the same site can interact. Since each of the M/2 identical 
sites can only support one bound state, the exchange integral 
must vanish and the Coulomb integral is the same for each 
interaction yielding £,-f/aLata/j,a;-|-. 

Since the Hubbard model and the Heisenberg model require 
magnetic fields to define problem instances, the electronic 
Hamiltonian which they consider is not precisely the same as 
equation ([6]); this also requires that the functional takes into 
account the separate spin components of the probability den- 



sity when performing the minimization in equation ( 17 1. 

Regardless, a polynomial time solution of UNIVERSAL 
Functional would also solve the QMA-hard electronic 
Hamiltonian with local magnetic fields. This follows as the 
functional is convex 



F lL Pi n i\ = m in Tl " 



(T e + V ee )pW 



^ LPJ ( min Tr (Te + Vee)p) 



Pj 



AN) 



(20) 



and the one-electron probability densities are also convex al- 
lowing one to perform convex optimization to find the min- 
imum energy of the QMA-hard electronic Hamiltonian with 



the local magnetic fields. Since the conditions for consistency 
are simple to check for one-electron probability densities^, 
the complications encountered in 2-RDM convex optimization 
from the approximate consistency conditions are not present. 

The inclusion of UNIVERSAL FUNCTIONAL in complexity 
class QMA follows along the same lines of the generic two- 
fermion problem discussed in section |3.4| However, in this 
case, Arthur must subtract the energy of the local field from 
the total energy to decide if the energy of F is above or below 
E T ±8. 

6.2 Restriction to pure states 

The pure state restriction of DFT affects the complexity of 
computing F. Consider evaluating F[n] where the input den- 
sity arises from optimization over pure states I^P) instead of 
mixed states as in equation ( fTTV In this case, the universal 
functional is no longer convex. 

£>, min (Wi\(T e +V ee )\%) ± min (V\(T e +V ee )\V). 

The optimization problem would then be NP and an oracle 
for Pure Universal Functional would allow a Turing 
reduction from QMA to NP. Stated differently, an NP machine 
with access to an oracle for Pure Universal Functional, 
could solve any problem in QMA. Similar to the 2-RDM case, 
the restriction to pure states complicates the verification such 
that two Merlins are required to show that the state is pure and 
that the proposed solution is correct. 

7 Other topics 

In this article, we have focused on three pillars of quantum 
chemistry but there are several other results worth mentioning. 
Brown et alM^ showed that calculating the density of states for 
quantum systems is in the same complexity class as computing 
the classical partition function. The complexity class of these 
problems is #P where the problem instances request the num- 
ber of solutions to problem in complexity class P. Rassolov 
and Garashchuk 20 briefly examined the complexity of the de 
Broglie-Bohm formulation of quantum mechanics finding that 
Monte Carlo techniques can be employed to efficiently sam- 
ple the propagated Hamiltonian-Jacobi equations without the 
quantum potential. They then suggest that the singularities of 
the quantum potential will ultimately be the source of difficul- 
ties in the quantum propagation. Troyer and Wiese 76 showed 
that the sign problem, properly defined, arising when evalu- 
ating quantum partition functions with Quantum Monte Carlo 
is NP-hard. Extension of DFT can be used in the time depen- 
dent domain in what is termed time-dependent density func- 
tional theory (TDDFT). Tempel and Aspuru-Guzik 77 recently 
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demonstrated that one can construct TDDFT functionals for 
use in quantum computation. 

8 Concluding remarks 

To close the paper, we make a few remarks on scaling of the 
error tolerances. In the computational problems investigated, 
the error tolerance was upper bounded by an inverse polyno- 
mial l/q(N) in the system size. In HARTREE-FOCK, it was 
the pseudo-spin mapping that necessitated the promise, but, 
in Universal Functional and N-Rep, the promise is in- 
herited from the 2-Local-Spin Hamiltonian problem. At 
first glance, it would appear that the error tolerance becomes 
more restrictive as the system size increases, however this is 
not necessarily the case. 

We focus on the problem ^-Local-Spin HAMILTO- 
NIAN since it is at the heart of the reductions for both UNI- 
VERSAL Functional and N-Rep. In this problem, the er- 
ror tolerance, 8, is normalized by the number of terms in the 
Hamiltoniarp 7 ] m. Therefore, the error tolerance only shrinks 
with system size if the number of non-zero terms in the Hamil- 
tonian is constant or grows slower than q(N). In many cases, 
the number of terms in the Hamiltonian increase with system 
size. Thus, even if asking for fixed error, e.g. 1 kcal/mol, as 
the system size increases the promise is still fulfilled and the 
problem of deciding instances of £-Local-Spin Hamilto- 
nian is QMA-complete. 

When the error tolerance is independent of the system size 
or grows slowly, it is not clear if the problem remains QMA- 
complete. This is related to the on-going research into possible 
quantum generalizations of the PCP theorem^] The PCP theo- 
rem shows that the task of approximating some NP-complete 
problems is also NP-complete'22122]. A quantum generaliza- 
tion, if it exists, would allow one to show that ^-Local-Spin 
Hamiltonian remains QMA-complete even when the nor- 
malized error tolerance is bounded by a constant instead of 
an inverse polynomial of N. This remains a prominent open 
question in modern computer science and is discussed more 
thoroughly by Osborne^. 

With respect to complexity in quantum chemistry, we have 
examined three computational complexity of different quan- 
tum chemistry algorithms: HARTREE-FOCK, N-Rep, and 
Universal Functional. The difficult of these computa- 
tional chemistry problems relative to problems found in other 



tt One may ask why not use the norm of the Hamiltonian? For a Hamiltonian, H, 
the operator norm is the maximum eigenvalue and computing the maximum 
eigenvalue is just as hard problem as computing the ground state of — H, thus 
also a difficult problem. 

%% The term PCP refers to probabilistically checkable proof systems where the 
verifier is only allowed to randomly access parts of a non-deterministically 
generated proof. Using only polynomial time and a fixed amount of random 
bits, the verifier is to validate or invalidate the proof. 



fields gives insights into why computational chemistry is diffi- 
cult. The difficulty of these problems imply that even a quan- 
tum computer will be unable to solve all instances since it 
is believed that BQP does not contain NP. However, even if 
quantum computers are never constructed, the study of quan- 
tum computational complexity gives new insight into the rel- 
ative difficulty of problems. For instance, characterizing valid 
two-electron reduced density matrices or evaluating the uni- 
versal functional of DFT is arguably more difficult than eval- 
uating the Hartree-Fock energy based on the probable sepa- 
ration of complexity classes QMA and NP. The examination 
of other properties of molecular systems and other problems 
encountered in quantum chemistry will provide fertile ground 
for future research into the difficulty of quantum chemistry. 
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